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Abstract 

We study 't Hooft's integral equation determining the meson masses M„ in multicolor 
QCD2. In this note we concentrate on developing an analytic method, and restrict our 
attention to the special case of quark masses mi = m2 = S'/a/tt. Among our results 
is systematic large-n expansion, and exact sum rules for M„. Although we explicitly 
discuss only the special case, the method applies to the general case of the quark 
masses, and we announce some preliminary results for mi = m2 (Eqs. (16.11) and (16. 3p ). 
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1 Introduction 



As was discovered by G. 't Hooft in 1974 [1], the mass spectrum of mesons in multi-color 
QCD in two dimensions admits for exact solution, because in this model the mesons are 
essentially the two-body constructs, and their masses are exactly determined by the Bethe- 
Salpeter equation. For the mesons built from two quarks of bare (lagrangian) masses mi 
and 1712, the Bethe-Salpeter equation reduces to the singular integral equation 



X 1 — X 



where 



«i = ^-l, a2 = —^-l, (1.2) 

g 9 

with g being the 't Hooft coupling constant (which in 2D has the dimension of mass). The 
function ^p{x) has to satisfy the boundary conditions 

y.(0) = y.(l) = 0, (1.3) 

whence Eq. fll.ll) defines the spectral problem for the parameter A; the eigenvalues A„,?t, = 
0, 1, 2, ... are discrete, and determine the meson masses. 

Ml = A„ . (1.4) 

In principle, the problem can be solved numerically, to any degree of accuracy, and over 
the years a number of approaches were developed to that end [1-5]. However, we believe 
equation (11.11) deserves further study from an analytical standpoint. In our opinion, the most 
interesting problem with respect to Eq. fll.ll) is understanding the analytic properties of the 
eigenvalues A„ as the functions of complex ai and a2- Without significant analytic input, 
straightforward numerical approaches seem to be unsuitable to addressing this problem. 
At the same time, the neat form of equation (11. ip suggests that perhaps some analytic 
information can be extracted. 

In this note we report new results about the spectrum {A„} in the special cas^ 

tti = ^2 = . (1.5) 
Among our results is systematic semiclassical (large-n) expansion of the eigenvalues. 



2A„ = n + ! - 



" 37r6 {n + ^ 



I 2(-l)"+^ [ 41og[vre^^-^(n + f)] ) .^^^ 



Here 7£; is the Euler constant and we display just three leading terms, but in principle any 
number of terms can be produced via our technique (next four terms can be deduced from 
(B^SD, fOTj) . fOHj) and equations (Q, ([13]) in Appendix A). Note the unusual logarithmic 
factors in the third and higher terms, which make this expansion look rather different from 



^Note that it is not the case of massless quarks. In particular, the chiral symmetry is broken. 
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the standard WKB expansion in the Schrodinger problem. In addition, our approach allows 
for analytic evaluation of the spectral sums 

oo ^ oo ^ 

with integer s = 2, 3, 4, . . . . (The sums here are over even or odd eigenvalues. The 
corresponding eigenstates are even or odd with respect to obvious x — > 1 — x symmetry of 
(11. ip .) For low s we have, explicitly 

= 7C(3), gL^) = 2, 

Gf = -|vr^ + 28C(3), gL'^ = -| + | tt^ , 

(1.8) 

Gf = -27r2 + 42C(3)-|7r^C(3) + f C^(3) + f C(5), 

G^^^ = ^-|^^ + ivr^C(3)-f C(5). 

Again, in principle analytic expressions for any given s can be obtained, but for larger s 
the calculations become increasingly involved. At the moment we have these numbers up to 
s = 13, but only those with s = 5, . . . 8 have sufficiently compact form to be presented in 
Appendix A. Put together, the large-n expansion (II. 6p and the sum rules (II. 7p provide good 
control over the entire spectrum: the large-m parts of the sums (ll.7p can be approximated 
by the asymptotic expansions (11.60 . thus providing equations for the lower eigenvalues. 

We regard this work as preparatory for studying the spectrum of (11.11) with arbitrary ai, 
a2, with the aim of understanding analytic properties of the eigenvalues at complex values 
of these parameters. We concentrate here on developing the technique, and the case (11.51) is 
convenient for testing its efficiency. Besides, many details have particularly neat form in this 
case. But for the most part, our technique admits more or less straightforward extension to 
the general case, which will be the next stage of this project. The method also seems to be 
suitable for analysis of a large class of Bethe-Salpeter equations of the type of (11.11) which 
emerge in many 2D field theories with confining interaction^. 

The paper is organized as follows. In Section 2 we discuss general properties of Eq. Ol.ip . 
In particular we relate its solutions to solutions of a certain functional equation (see Eq. (l2.6p ) 
of the type of Baxter's T — Q equation, with special analyticity. In Section 3 we develop 
A-series expansion of the solutions of this equation. This expansion generates analytic ex- 
pressions for the spectral sums (11.70 . Asymptotic expansion at A ^ cxd is developed in 
Section 4. It results in the large-n expansion of the eigenvalues A„. In Section 5 we test 
these results against the numerical solution of (11.10 . 

While this paper was in preparation, we have made some progress in studying the more 
general case of (11.10 . with nonzero but equal values of the parameters 

ai = 02 = a . (1.9) 

^This situation is typical when one takes a field theory with exact vacuum degeneracy and adds small 
interaction which lifts the degeneracy, giving rise to the confining force between the kinks. The simplest 
example is the Ising field theory, in the low-temperature regime, in the presence of a weak magnetic field [6,7]. 
Unlike the multicolor QCD, where the equation (|l.ip is exact in the limit Nc ^ oo, in that case the associated 
Bethe-Salpeter equation is only an approximation, expected to be valid when the magnetic field is sufficiently 
small, but it seems to produce meaningful insight into the mass spectrum even at a large magnetic field. 
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We intend to devote a separate paper to discussing this more general case, where indeed a 
very interesting analytic structure of A„(a) emerges. But we could not resist the temptation 
to announce some results here, which are presented in Section 6. 

2 Functional equation 

We find it useful to recast Eg. (11. II) into somewhat different form, via the integral transfor- 
mation 

f"^ (\v dr 
V(^)= ^^(^)(t^)^' ^(^)= / . Vix) fe)"^, (2.1) 



2n ' ' ^i-^'^ ' Jo 2x(l -X 

which is just Fourier transform with respect to the variable | log{j^) (This transformation 
was previously used in Ref. [8]). The u-space form of (11.11) is 

/oo 
dz/' ^(z/-z/') ^(z/') = 0, (2.2) 
-oo 

where the kernel 

in the right hand side is regular at all real u. The solution \E'(z/) must decay at |z/| oo 
(for the norm || (p ||^= dx \(p{x)\'^ to be finite), and it must be a smooth function of u (for 
the function ip{x) in (11. ip to satisfy the boundary conditions (11. 3p ). The spectrum {A„} is 
determined by the existence of solutions which satisfy these conditions. In fact, both these 
conditions, once satisfied, are satisfied with substantial redundancy. 

Equation (12.21) dictates that any smooth solution is in fact analytic. Moreover, it is 
possible to show that the solutions ^E'(z/) are meromorphic functions of z/, with the poles at 
u = ±{2k — 1) i. A; G Z, of the order A; G N. In particular, the function Qlu) defined as 

Q(z/) = z/ cosh (f) ^(z/) (2.4) 

is analytic in the strip |53mz/| < 2, grows slower then any exponential of z/ at infinity, and 
turns to zero at z^ = 0, ±2i, i.e. 

Q(0) =Q(±2i) = 0. (2.5) 

Under these conditions the integral operator in the right hand side can be inverted in terms 
of finite difference operator (which is derived by standard manipulations with shifts of the 
integration contour), leading to the functional equation 

g(z/ + 2i) + g(z/ - 2i) - 2 g(z/) = -47rA z/^^ tanh (^) Q{u) . (2.6) 

Equation (12. 6p is the basis of our analysis below. 

A quick look at the asymptotic form of (12. 6p at 3?e z/ — >■ cxo reveals that its solutions gen- 
erally behave as e^" f{v), with integer /c, and /(z/) bounded by any exponential. Obviously, 
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any positive k would violate the asymptotic condition for \I'(z/). Thus, we are interested in 
the solutions which are bounded as 



Q(zy) = 0(e^l'^') as |SRez/|-.oo, (2.7) 

with any e. Note that this condition implies that the function "^{ly) in fact decays exponen- 
tially in this limit. 

A solution of ( 12. 6p with the desired analytic and asymptotic properties exists only at 
specific values of A, which determine the eigenvalues of (12. 2p . However, if the conditions (12.51) 
are relaxed, the solutions Q(z/|A) exist at any A. For generic A, the associated function \l/(i/|A) 
no longer satisfies the integral equation (12.21) . Instead, it solves the related inhomogeneous 
equation 

/oo 
di^' Sii^-i^')<iJ{u'\\)=F{i^\\), (2.8) 

where 



J — c 



with the coefficients ^^-(A) and q-{\) related to the values Q(0|A) and Q{±2i\\) in a linear 
manner (note that in view of the functional equation (12. 6p . only two of these values are 
independent). Since now in general Q(0|A) ^ 0, the integrand in the l.h.s. involves a first 
order pole at z/' = 0, and the integral is understood as its principal value. It is possible to 
show that, given the coefficients q±, the solution of (12.81) is unique. These coefficients can be 
chosen at will, and therefore the equation (12. 8p generates two-dimensional space of functions 
\l/(z/|A). It is natural to choose the basis in accord with the obvious u — —u symmetry of 
the problem. We thus define symmetric and antisymmetric basic functions, 

vl/±(-z/|A) = ±vl/±(z/|A), (2.10) 

which solve the equation (12.81) with F(z/|A) in the r.h.s. taken to be 



--^v-y w^,.x and F^(u) = — , (2.11) 
' sinh(^) ^ ' sinh(f ) ' ^ ' 

respectively. 

At the spectral points A = A„ the original equation (12. 2p is to be recovered. That means 
that at certain values of A the basic functions \I'±(z/|A) diverge. More precisely, Eq. (12. 8 p can 
be rewritten in the form of an inhomogeneous Fredholm integral equation of the second kind 
(see Appendix B for details) and it follows from the resolvent formalism that \l/-|-(z/|A) are 
meromorphic functions of A, with only poles at the eigenvalues of (12. 2p . 



A — A2m A — A2m+1 



m=0 m=0 
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where, as we have mentioned in Introduction, and A2m+i, m = 0, 1, 2, . . ., refer to the 
eigenvalues of (12.21) in the even and odd sectors, respectively, and '^2m{^) and '^2m+i{^) are 
associated eigenfunction^. 

It is useful to note that the functions \E'+(z/|A) and ^_(z/|A) are related to the "quark 
form factors" of the vector current = i^'Ji^ip and the scalar density S = ipip, respectively, 
with the parameter A (more precisely 2tt g'^ X) interpreted as g^, the square of the total 2- 
momentum (see Refs. [9,10], where inhomogeneous integral equations equivalent to (I2.8p . 
(12.111) appear in this connection). Therefore the structure (12.121) is well expected, and the 
coefficients c„ in (12.12^ are related to the matrix elements 

(0 I J^(0) I M2^,g) = ie^^g'^v^vri C2„ , (2.13) 



(0 I 5(0) I M2m+i,g) = 2T!:g a/A^ C2m+i , 

where | Mn,q) stands for the n-th meson state with 2-momentum q. Let us mention here 
the neat expressions for the current-current correlation function in terms of \I'+(z/|A), 

{Mq)M-q)) = '-^ {^-9,.) [l-^+m] ■ (2.14) 

Having in mind this analyticity in A, our strategy in solving the problem will be as follows. 
Starting with equation (12.61) . we will be looking for two solutions, Q+iylX) and Q-{u\\) of 
the functional equation (12. 6p . analytic in the strip |53mz/| < 2, and growing slower than any 
exponential at |3fJet/| oo. We will assume that 

Q±{-v\\) = tQ±{i^\X), (2.15) 
and fix the normalizations by the conditions 

g+(2i|A) = -g+(-2i|A) = 2i , g-(o|A) = i . (2.16) 

Then the functions \I'-|-(z/|A) related to g±(z^|A) as in (12. 4p have appropriate symmetry (I2.10p . 
and solve the equation ( 12. Sp precisely with the right-hand sides ( 12. lip , as one can readily 
verify. In fact, the remarkably simple formula 

dx\ogD±{\) = 2id, \ogQ^{v\X)l^. (2.17) 

relates the logarithmic derivatives of at = i to suitably defined spectral determinants, 

m=0 ^ ' m=0 ^ ' 

(To be precise, somewhat more complicated expressions, Eqs. (18. 20 p . follow directly from the 
integral equation (12. 8p . See Appendix B, where we explain the status of Eq. (12.170 .) Note 
that this relation is insensitive to normalization conditions assumed for Q±[y\\)- In what 
follows, we develop two different expansions for such solutions Q±{y\S)- One is just the 



■^Here and below ^'„(i^) stand for normalized eigenfunctions, i.e. we assume that jj^ da; |iy9„(a;)p ~ 1 for 
the associated 1^5^(2;), Eq. (|2.ip . 
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power series in A, and the other is the asymptotic expansion around the essential singularity 
X = oo. Eq. fl2.17p translates the former expansions into the sum rules i \1.7\i . while the latter 
ones lead to the large-n expansion (11.61) . 

Before turning to the details, let us make the following remark. The functional equation 
(12. 6p has the form of the famous T ~ Q relation of Baxter, and many general statements 
can be adopted to our case. In particular, it is easy to show that the so-called quantum 
Wronskian built from the two solutions Q±{i'\X) is a constant, 

Q+(z/ + i|A) Q4u - i|A) - g+(z/ - i|A) g_(z/ + i|A) = 2i . (2.19) 

The fact that this combination does not depend on u follows from the functional equation 
(12. 6p . and the asymptotic conditions at |3?ez/| — ^ oo. The particular value 2i in the r.h.s. 
reflects the special normalization (12. 16^ : with different normalization it would be a differ- 
ent (generally A-dependent) constant. This equation, combined with (12.171) . allows one to 
establish some useful relations. It follows from (I2.17P that (5+(i|A) turns to zero at the odd 
spectral values A = X2m+i (likewise, Q_(i|A) does the same at the even values X2m)- But the 
identity Q+(i|A) (5-(i|A) = i (elementary consequence of (I2.19p ) shows that X2m+i exhaust 
all zeros of (5+(i|A) viewed as the function of A. In other words, (5+(i|A) must be proportional 
to D_(A)/Z^+(A), up to a factor which is entire function of A with no zeros, i.e. in our case 
the factor of the form exp(a + bX). More careful analysis (in the next section) allows one to 
fix this ambiguity completely. Let us present here the result in the form 

Q+(i|A) ^ 1 lUXl Q4i\X) ^ ZMA) . . 

g+(2i|A) 2 D+{x) ' g_(o|A) D_(A) ' ^ ■ ' 

insensitive to normalizations of Q±{i^\X). 



3 Expansion in powers of A 

In principle, one can generate the expansion in A just by iterating the integral equation (12. Sp . 
with the right-hand side taken in one of the forms (12. lip (see Eqs. (I8.14p in Appendix B). 
This leads to the convergent series 

oo 

Q^{u\X) = J2Q±\^) (3-1) 

s=0 

with the coefficients given by s-fold integrals involving the kernel 5'(z/). Direct evaluation of 
these integrals is difficult, and therefore we take another approach based on the functional 
equation (12. 6p . We look for the solution of (12. 6p in the form of the power series (13. ip . with 
the coefficients Q^^\i') having the symmetry (I2.15p . analytic in the strip |S5mz/| < 2, and 
growing slower than any exponential at |3fJez/| —>■ oo. It is clear upfront that at each order 
these conditions fix the coefficients uniquely (after all, it is just a somewhat indirect way of 
evaluating the integrals appearing in the iterative solution of (12. 8p ). But the solution for 
Q±\^) obtained this way involves polynomials of u of growing degree, and the expressions 
quickly become cumbersome. The following observation greatly facilitates the calculations. 
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Note that the factor 



27rA tanh (^) 



(3.2) 



in the r.h.s. of this equation, is insensitive to the shifts z/ ^ z/±2i, and if no attention to the 
analytic properties is paid, it can be regarded as constant. Then the equation (12. 6p . written 

as 



Q{u + 2i) + Q{u - 2i) - 2 Q{u) = -2z u'^ Q{u) 



(3.3) 



is recognizable as one of the recursion relations satisfied by confluent hypergeometric func- 
tions [13]. Specifically, the functions z/M(l + f , 2, -iz) and 1(1 + f ) f/(l + | , 2, -i^) are 
known to satisfy (13. 3p (see e.g. [13], Eqs. 13.4.1, 13.4.15). Here the conventional notations 
M(a, c, x) and U (a, c, x) for two canonical solutions of confluent hypergeometric equation 
are used. Of course, by themselves these functions do not provide solution to our problem, 
since they have wrong analyticity in v. For one, the second of these functions has logarith- 
mic singularity at 2; = 0, which in view of (13. 2p produces unpleasant branching points in 
the i/-plane. This problem is easy to cure by observing that the logarithmic term by itself 
satisfies Eq. (13.31) . and subtracting it produces another solution which is now a single- valued 
function of u. Thus, we found it convenient to use the combinations 



M.{u,z) 



z/e^ M(l + f, 2,-iz) , 

-ize^r(l + f)[/(l + f,2,-i^) - 
l[z\og{-\ze'^^)+i7i^X] M+{u,z) 



(3.4) 



(3.5) 



where the coefficients e'2 and the extra constant in the brackets in (13. 5p are chosen to ensure 
the symmetry 



M+(-z/, -z) 



TM±(z/, z) 



(3.6) 



in accord with the obvious symmetry of Eq. (13.30 . Both (13. 4p and (13. 5p are entire functions 
of z, in particular both can be represented by convergent expansions in the powers of z: 



where 



Mju^z) 



f ^ (1 + f ) 



s=Q 



(s + l)\s\ 



(3.7) 



M_(z/,z) = i [e^S(i^,z) +e"^S(-z/, -2) 



E(z/, z) 



^ 2 ^ 



s=l 



s! s-1 ! 



^(s + I) - V;(s) - ^Pis + 1) + ^(i) (-i^)^ . (3.8) 



These expansions make explicit a more serious problem. In view of (13. 2p . each term of this 
expansion produces poles at z/ = ± i, of growing order, and thus both ( 13.40 and ( 13. 5p . viewed 
as the functions of u at fixed A, have essential singularities at these points, whereas we need 
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solutions of (12. 6p analytic in the strip |53mi/| < 2. We are thus compelled to look for the 
solutions in the form 

Q±{iy\X)=A±{z,X) M±{u,z) + B±{z,X) zM^{u,z) , (3.9) 

where the coefficients, entire functions of z^, are to be adjusted to compensate for the above 
singularity at z = oo. So far we were unable to find the closed form solution of this analytic 
problem. But it is easy to generate the solution as an expansion in the powers of A. In view 
of the above analyticity, we assume that the coefficients can be expanded in double series 
in A and z"^. Regarding them as the functions of u and A (through the relation (13.21) ). this 
expansions has the form of the power series 

oo oo 

A^{z, X) = J2 4\r) X^ , B^{z, X) = J2 b^±ir) (3-10) 

s=0 s=0 

with a^±{T) and &±^(t) being polynomials in r = (^)^ = x tanh^(^), of the highest degree 
[s/2]. The numerical coefficients in these polynomials are to be adjusted in such a way as to 
compensate for all the pole terms generated by the expansions of the functions M±{i>, z) in 
(13. 9p . order by order in A. The remaining constant terms are then fixed by the normalization 
conditions (12.161) which demand that ^±(0, A) = 1. Clearly, this linear problem at each order 
has a unique solution. We have calculated explicitly the polynomials 

«±V), b^±ir) up to 
s = 13. Let us present here the first few of them, just to give the fiavor of it: 

J^)-^r, aj) = ir ( 50 + 21 C(3) - 5r ) , 
6« = |, 6(2) = 1 (6 + 7C(3)-r), (3.11) 



— T , a_|_ — g 

= i , ^^^^ - ^ ^^^^ 



and 



aL'^ = -r, aL'^ = |r, aL^^ = ^ r ( 21 C(3) - 14 + 3r ) , 

&L°^ = 0, b^l^=4, 6^'^ = IC(3) -5-r. (3.12) 

Eq. (l2.17p makes it straightforward to convert the A-expansions of Q±(z^|A) into the ex- 
pansions of the spectral determinants (I2.18p . 



logD±(A) = {lTl)X-Y,s-' A^ , (3.13) 

s=2 



where the coefficients give explicitly the spectral sums (II. 7p . For the few lowest s the result 
of this calculation was already displayed in (II. 8p . but we present many more in Appendix 
A. With many Ci^^ known, the sum rules ( 11. 7p become a useful tool in determining the 
eigenvalues A^, especially so when combined with the large-n asymptotic expansions, which 
we derive in the next section. 
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4 Asymptotic expansion at A ^ oo 

To develop the large-A expansions of the functions Q±{i^\\) we start by constructing a formal 
solution of the functional equation (12.61) . of the following structure 

oo 

5(HA) = (-A)-T J2s,{u)\-'. (4.1) 

It is impossible to satisfy all the analytic conditions required for the functions Q±{iy\X) 
within this ansatz, but we would like to get as close to the desired analyticity as possible. In 
particular, we demand that the coefficients S'fc(z/) are meromorphic functions of z/, growing 
slower then any exponential at |3fJez/| oo. The form (14. ip is obviously designed to serve 
the case of negative real A, if we choose the principal branch of (— A)^'^ (other branches 
exhibit unacceptable exponential growth at |3fJez/| oo). 

Plugging this expansion into (12. 6p generates a sequence of recurrent functional equations 
for the coefficient functions ^^(z/). In the zeros order we have 

So{u + 2i) = ^ tanh (^) So{iy) . (4.2) 

The solution of this equation, analytic in the strip |53mi^| < 2 and bounded at |3ftez/| oo, 
is unique up to a normalization. It can be written in an explicit forrr0 

in terms of the Barnes G-function (see e.g. [14]) 

oo r- 

G(x+ 1) = (27r)2 e ^ [[ y- + - 



n=l 



e-"+f^^ 



(4.4) 



At higher orders in A ^ the equation (12. 6p leads to the recurrent relations of the form 

afc (z/ + 2i) - (Tfc(z/) = pfc(z/) (4.5) 



*The function ipoiv) — ^ tanh(^) Sq{v), with 5*0(1^) as in (|4.3p . provides an exact solution to the 
"scattering" problem 



associated with (jl.ip . see Ref. [11]. Namely, 

(p{x) = / dt^ ^Jq{i^) ■ 



Using the known asymptotic behavior of the Barnes G- function [14], it is straightforward to derive the 
"scattering phase" in 

ip{x) c^'"^ + e^^ e'"^ as a: ^ 00 . 

from which the constant term | in Ea. (|1.6p (already conjectured in Ref. [1]) follows. Our analysis in this 
section goes beyond this simple approximation. 
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for the ratios a"fc(z/) = Sk{u) / Sq{u), with Pkiy) being certain expressions involving 0-^1 {u) 
from the lower orders k' = k — 1, k — 2. While beyond the leading order no solutions analytic 
in the strip |53mz/| < 2 exist, it is possible to find solutions analytic in that strip except for 
the points z/ = 0, ±2i, where the poles of growing order appear. The result of this calculation 
is summarized by the formula 

S{u\X) = R{z,X)U{u,z) , (4.6) 
where U{i', z) stands for the formal asymptotic series 

U{u^z) = (-A)-- Soiu) f2 ^^^^^y^^' (4.7) 

fc=0 

and z is the same combination (13. 2p . insensitive to the shifts u ^ u ±2i. The fact that (14.71) 
satisfies (12. 6p can be verified directly, but it is clear upfront from the following observation; 
The series appearing in (14.71) . when multiplied by r(l + ^) (—iz)~^~''^ , coincides with the 
asymptotic expansion of the function r(l + y) U{1 + ^-^,2,—iz), which satisfies (12. 6p . In 

writing (14.71) we simply replaced the overall factor r(l + y) {—iz)^^^'~^ by the much more 

analytically attractive (— A)^^ 5'o(z^). The factor R{z,\) in (14. 6 p represents the ambiguities 
in the solutions of (14.51) : it is to be understood as a formal series in the powers of z~^ and 
A~^, or equivalently as a series in A~^ with the coefficients being polynomials in the variable 

c = ivr coth (^) . (4.8) 

The asymptotic expansions of the functions Q±{i'\X) can be built from the formal solution 
(14. 7p in much the same way as the A-expansions were constructed from the basic functions 
(13. 4p . (13. 5p in the previous section. Having in mind the symmetry (I2.15p . we look for Q±{'y\X) 
in the form[fl 

Q±{iy\X) X R±{z, A) U{u\z) t R±{-z, A) U{-u\ - z) . (4.9) 

The coefficients i?±(z. A) are to be adjusted to fix the analytic problems present in lJ{y\z) 
and lJ(^—v\ — z). One of these problems was already mentioned above. The series (14.71) 
explicitly exhibits at each order in A~^ poles at u = 0, ±2i, of the growing order. This 
problem can be fixed order by order in A^^, with R±{z, A) taken in the form 

00 

R±{z, X) ocl + Y, R±\c,L)X-'' , (4.10) 

k=l 

with R^\c,L), polynomials in the cotangent (14.81) . adjusted to cancel these poles. Because 
of the factor (— A)^~, the Laurent expansions of (14.71) around the poles generate logarithms 
of —A. As a result, the coefficients R^^^ in (14.101) emerging in this calculation turn out to be 
also polynomials in the variable 

L = log (-27rA) + 7s . (4.11) 
^Here and below the symbol x stands for equality in the sense of asymptotic series. 
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This is a novel feature of the large-A expansion, which ultimately leads to the logarithmic 
factors in the expansions (11.61) . As expected, the solution of this pole-cancellation problem 
at each order in turns out to be essentially unique, i.e. unique up to terms which 
can be absorbed into the overall normalization of Q±{i'\X). With the relations ( ]2.20p . the 
normalization conditions (12.161) imply the following general form of the coefficients R±{z, A): 



-A) 2 i^ 





(A) 1 


±1 r 


2D. 


f(A)J 





(c,L) A- 



k=l 



(4.12) 



We have explicitly computed the polynomials P^\c,L) up to A; = 7, but display here only 
the first few of them (again, just to give the fiavor of the emerging expressions): 



Pi'\c,L) = 0, P^'>{c,L) 



(2), 



±1 

4^ 



(c,L) 



±1 
24^ 



(6c- 12L + 6T 1) . (4.13) 



The overall factors in (I4.12p are related to the ratio of the spectral determinants (I2.18P via 
(12.201) : the expansion 



exp 



L 



2A 



2(A7r 



2^2 



12(A7r 



2^3 



(4.14) 



is obtained in a straightforward way once P^'' (c, L) are determined, by imposing the nor- 
malization condition (I2.16P order by order in A^^. 

These results are readily applied, through Eq. (l2.17p . to write down the large-A expansions 
of the individual spectral determinants 



9AlogD±(A)xL-l + log(4) + 



-1 ±2 
8A 



TT 



5^pf (0,L) A- 



(4.15) 



k=2 



where L is the logarithm (14.110 . and the terms oc A~^ and higher involve the polynomials 
(I4.13P specified to c = 0. This equation determines the large-A expansions of D±{X) up to 
an overall numerical factors 



D±(A) xrf± (87re-2+^^)^ (-A) 



8 4 exp 



5^pf(L)A- 



k=l 



(4.16) 



where the polynomials F^\l) are easily deducible from (I4.15p . e.g. 



47^2 



1 ± 12L 

487r4 



etc 



(4.17) 



One immediate consequence of the asymptotic expansions (14.160 are analytical predictions 
for the regularized sum 



G 



(1) 



m=0 



1 

A2m 



1 

m + 1 



1 



1 



m=0 



A2m+1 m + 1 



(4.18) 
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The form of the pre-exponential factor in fl4.16l) imphes 

G'i^ = log(87r) - 1 , G'l^ = log(87r) - 3 . (4.19) 

The numerical factors d± can not be obtained from (14.151) . In fact, at the moment we do 
not have analytic expressions for these constants. However, the exact relation 

follows from (14.141) . Note that the constants d± can be written as (fast convergent) products 

' m=0 ' m=U 

The relation (I4.20p is a rather nontrivial prediction of our theory. Fortunately, the constants 
d± play no role in derivation of the large-A expansion of the spectrum below. 

It is important that the form (14. 9 p was designed to describe the asymptotic behavior 
of Q±{i^\X) at large negative A, therefore (I4.15P generates the asymptotic expansions of the 
spectral determinants at A — >■ — oo. In view of the analytic structure (I2.18p . these expansions 
are in fact valid at all (sufficiently large) complex A, except for when A lies in a narrow sector 
around the positive real axis in the complex A-plane. But since the main object of our interest 
is the spectrum {A^}, we are especially interested in the asymptotics of D±(A) at real positive 
A. Below we argue that the asymptotic behavior in this domain is correctly described as 

D±(A)xZ}W(A) + d1-)(A), (4.22) 

where D^\X) and \X) are the results of term-by-term analytic continuations of the 
series (I4.16P from the negative to the positive part of the real axis, in the clockwise and the 
counterclockwise directions, respectively (formally, D^\X) = D±{—e'^'^X) and D^\X) = 
D±{—e^'^X), where D±{X) are understood as the series (14.160 ). Then we can write the A 
+00 expansions as 

D±(A) x2rf± (87re-2+^^)^ A^-5±i e"±(^) cos [ f ( 2A - i ± | - $±(A) ) 1 , (4.23) 
where S-i-(A) and $±(A) are the asymptotic series of the form 

oo cx> 

s±(A) = Y.E'i^i) A-^ $±(A) = Y,<^t\i) a-^ (4.24) 

k=l k=2 

Here the coefficients and ^±\l) are polynomials in the real logarithm 

/ = log(27rA) + , (4.25) 
directly related to the polynomials F^\l) in (I4.16p . 

S^)(0 = i [Ff (/ + ivr) + Ff (/-iTT)], (4.26) 



$J)(0 = 1 [Ff (/ + ivr)-Ff (/-ivr)]. 

TT 
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Of these, $±(0 especially important since they enter the "quantization conditions" 

2m, (4.27) 

2m + l, (4.28) 

k=2 

which, with m = 0, 1, 2 . . . , determine the eigenvalues and X2m+ii respectively. There- 
fore we present explicitly $±^(/) up to = 7 in Appendix A (see Eqs. (17.31) and (I7.4p ). The 
large-n expansion (11.61) follows directly from (I4.27p . (14.281) . 

At the moment we do not have completely satisfactory proof of (14.221) . However there is 
a body of supporting arguments. The most important concerns the behavior of the functions 
Q±(z/|A) themselves at large positive A. The easiest way to understand the situation is again 
through the analytic continuation in A. The expression (14. 9 p can be analytically continued 
to positive A term by term in the expansions (14. 7p and (14.120 . With any such continuation, 
it still satisfies the functional equation (12.60 order by order in A, and its coefficients are still 
free of poles in the strip |53mz/| < 2. But there are two natural ways of the continuation - 
one is through the upper half-plane, and another is through the lower one. Thus at positive 
A we have two series-like solutions of ( 12. 6p . with correct analyticity in the strip IQ'mz/l < 2, 
both for (5+ and Let us denote them Q''^\v\\) and ''('^|A). The problem is that 
each of them exhibits unacceptably rapid growth at |3?ez/| ^ 1. It is possible to show that 
at, say, positive z/ — > +oo they behave as 

Q^+\u\\) ^ e-("+i) i?±(-27rA,A)|i=,_i^ e^^ M(z/, 27rA) , 

Q^±\u\\) ^ ±e-^-(^+i) /2±(27rA,A)U=i+i^ e"'^ M(z/, 27rA) . (4.29) 

Here M is a certain combination of the hypergeometric functions M±{v, 27rA) (13. 4p with co- 
efficients which, unlike exponential factors e^''^^, admit large-A expansion similar to (I4.10p . 
Note that this behavior is completely compatible with the functional equation, but contra- 
dicts the required large-i/ behavior of true functions Q±{v\X)i which must grow slower then 
any exponential. Admittedly, we are dealing here with asymptotic series in A~^, and using 
them to judge the v ^ oo asymptotics is problematic. But the series in (14. 7p is expected 
to be approximative at large A as long as \/~\ ^ v (and even more so for the series (I4.10p ). 
If one focuses on the region -\/A ^ z/ ^ 1, the exponential growth (I4.29P is clearly incom- 
patible with the expected behavior of true functions Q±(i^|A), which at z/ 1 must quickly 
(with exponential accuracy) become linear combinations of the functions M_|_(z/, 27rA) and 
M_(z/, 2'7rA) defined in (13.41) . However, it is clear from (I4.29P that one can form special linear 
combinations of Q^^^ and Q*-"-* in which the growing terms cancel out. These combinations 
involve the factors e^'^'^ and e~"^^ which do not admit A~^ expansions; this is why straightfor- 
ward A~^ expansions are impossible at positive A. The most compact way to describe these 
linear combinations is in terms of somewhat differently normalized functions 

Q±(z/|A) = 2±i i^ D±(A) Q±(z/|A) ; (4.30) 

instead of (I2.16P they satisfy 

Q+(2i|A) = 2 D+(A) , Q_(0|A) = i D.{\) . (4.31) 



k=2 

oo 

2A-|-5^$L'^(/) A-'= = 
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The normalization factors in fl4.30p make Q±{jy\X) entire functions of A. In particular, instead 
of (12.121) . at the spectral values of A we simply have 



'^2mM oc 



2m 



V cosh( 



2 > 



Q-{v\\ 



V cosh( 



2 > 



(4.32) 



At negative A (indeed, at all complex A except for the narrow sector around the positive 
real axis), the large-A expansion can still be written in the form (14. 9p . with the coefficients 
R±{z^ A) replaced by 



7^±(;z, A) = (8 



vre 



l + 5^7^^)(c,L)A 



(4.33) 



with new polynomials 7l±{c,L) which are obtained by combining (I4.12p with (I4.16p . Now, 
let Q±~\i^\)^) and Q^. \i'\X) be two asymptotic expansions obtained by formal term-by-term 
analytic continuation in A from negative to positive A, one through the upper half-plane and 
another through the lower one. It turns out that it is exactly the sums 

Q±^ + Q± ^ in which 

the unacceptable growing terms (14.291) cancel out. Thus it is natural to assume that correct 
asymptotic behavior of true functions Q±(z^|A) at real positive A is given by these sums[§. 



Q±(^|A) 



A — s> +00 , 



(4.34) 



Prom this, the form (14.221) immediately follows. We note here that after the cancellation of 
the growing terms, the combinations (14.341) have the following behavior at real |z/| ^ 1, 



Q+(z/, A) ~ M+(|z/|, 27rA) , Q_(i/, A) ~ sgn(i/) M+(|z/|, 27rA) . 



(4.35) 



It turns out that these equations give very good approximations of the functions even at 
1/ ~ 1, and even if A is not particularly large (see next section). 

Another piece of evidence supporting (I4.22p is numerical. First of all, the numerical values 
of X2m and X2m+i obtained from (14.271) and (14.281) . with some reasonable number of terms 
in the A~^ expansions included, provide remarkably accurate estimates for the eigenvalues, 
even for low levels. We discuss these numerics in the next section. But one can also match 
the large-A expansions of the spectral determinants to the power series expansions 



D±(A) = l + (lTl)A + ^ Di^) A^ 



(4.36) 



The latter converge in the whole A plane. With many terms included, the expansions (14.360 
are expected to approximate the functions -D±(A) well even if |A| is not small. Since we 
know as many as 13 terms of ( 14. 360 . one expects to have substantial domain at negative 
A where the truncated series (I4.36P match the asymptotic expansions (I4.16P (again, with a 
reasonable number of terms in the sum). More crucially, if (14.221) is correct, there must be 
substantial domain of positive A where it matches (I4.22p . This comparison requires knowing 



^The situation is reminiscent to how the WKB expansions of the wave- functions in quantum mechanics 
are matched around the turning points. 
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Figure 1: Plots of small- and large- A expansions of e~^'^^ The A-expansion, with terms 
up to oc A^^ in fl4.36p . is shown as the dashed line. The solid lines represent the large- A 
expansions, i.e. fl4.16p at negative A, and fl4.23p at positive A; in both cases terms up to 
oc A~^ are included. 

^-2.5 A ^_ 




A 



Figure 2: Same as in Fig. 1, but for e '^■^'^D^. In this case the small-A expansions is 
truncated to terms oc A^^. 

the constants d± in Eqs. fl4.16p . (14.220 . We use here the numerical estimates from the next 
section (see Eq. (15. ip ). In FigsiTl[2]we present simultaneous plots of the A expansions (I4.36p . 
with as many terms as are available, and the large-A expansions (I4.16P and (I4.23p . with 
the sums including all terms up to oc A"^. In fact, the plots are for e~^'^^ i5±(A), with 
the exponential factor added to make interesting parts of all three plots visible in the same 
picture {D±{\) themselves develop large amplitudes already at A ~ 1). As expected, there 
is a good match at negative A between —1 and —0.4, but one can also see a clear match 
at positive A, in the domain between 0.6 and 1.6 where the functions already show "live" 
behavior. Note that the two lowest zeros of both Dj^ and -D_ are already visible at these 
orders of the A-expansions. In fact, the positions of the lowest zeros Aq and Ai stabilize 
rather fast as one adds more and more terms to (I4.36p . This convergence is particularly 
impressive for Aq. Fade approximation of the A-expansion of the ratio ^^4^, Eq. (l2.20p . 
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yields the following estimate of the lowest eigenvalue, 

2Ao = 0.737061746292690 
Compare this number to the numerical result in Table 1. 



(4.37) 



5 Numerical results 

As was mentioned, it is not difficult to compute eigenvalues {A„} by direct numerical solution 
of equation (11. ip . A variety of numerical methods exists [1-5]. We have used the expansion of 
fix) in Chebyshev polynomials from [2,4], which seems particularly suitable in the case ( 11. 5p 
since it automatically guarantees the function (p{x) correct behavior near the boundaries x = 
0, 1 (besides, its implementation requires perhaps the least amount of creative programming). 
With this method, fourteen significant digits for as many as 50 lowest eigenvalues A„ can be 
obtained by truncating to matrices of the size 400 x 400. Below we use the notation An™™'' 
for these numerical estimates. 

In Tables 1 and 2 we compare these numbers, for even and odd n separately, with the 
results of large-A expansions. The first column in each of these tables shows numerical 
values yielded by Eq. fll.6l) . with all terms explicitly written there included. One can observe 
significant improvement as compared to the leading semiclassical approximation A„ ~ n + |, 
even for the low levels such as Ai and A2. The approximation (II. 6p corresponds to truncating 
the sums in Eq. (l4.27i) and (14.281) to terms ~ A~^, but one can obtain further corrections by 
including higher-order terms. We denote X^^ the estimates from equations (14.271) . (14.281) 
with terms up to oc A^'^ included, and present the numerical values of An''"' (together with 

/iy\ /iy\ 

the deviations 5A„ = An — An to show the expected accuracy of this approximation). 
Since we are dealing with an asymptotic expansion, one does not expect it to work well for 
low levels, but Tables 1,2 show that including these further corrections results in noticeable 
improvement even for levels as low as A3 and A4, and for higher levels the improvement 
becomes impressive. For n > 30 X^J^ are indistinguishable from An""™"* within the accuracy 
of the latter. 

Another impressive agreement is in terms of the sum rules (ll.7p , (I4.18P . One can evaluate 
the spectral sums in (11.71) . (I4.18P using the numerical values An""™"*, and compare these nu- 
merical estimates [ (^|^) j -with the analytic predictions (II. 8p . (I4.19P and (17. ip . (17. 2p . In 
fact, for low s the sums do not converge that fast. For instance, to estimate [G^f^ j(num) 
fourteen digits one needs to include as many as 10'' eigenvalues. Of course, this problem is 
easy to solve since we have very good large-n asymptotic approximations. In the sums (II. 7p . 

starting from some sufficiently large n one simply replaces A^T"™^ by the asymptotic form, 
(7) 

say An • In Table 3 we show numerical estimates obtained in this way for s = 1, 2, . . . , 8. 
It is an easy and pleasant exercise to check that these numbers agree with the analytic ex- 
pressions (II. 8p . (I4.19P and (17. ip . (17. 2p to all digits presented. As was mentioned, we actually 
have analytic expressions for Gj!^ with s up to 13, and we have verified similar agreement 
for these higher values of s as well. We also have computed the products (I4.2ip with the 
numerical spectrum, 

^(num) ^ 0.963178456398 , dL"''"'^ = 0.433582639833 . (5.1) 
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n 


zAr, irom ii^q. (11. ol) 




2 5Ai^^ 


r, (num) 
ZAn 


U 


U.73U 








0.73706174629269 


2 


2.748145 


2.748159 


6.3 


X 


10^ 


-5 


o '7/ioi/^nnioo '-7f\n 

2.7481609123706 


4 


4.749299 


4.7492955 


1.8 


X 


10" 


-6 


A '^/inonrooino^r 

4.7492953810375 


6 


6.749631 


'7 /I r\/^or\ /I o 

6.74962943 


1.7 


X 


10" 


-7 


6.7496294196488 


o 

8 


8.7497729 


O '7/ir\'7'71 r O /I 

8.749771584 


2.9 


X 


10" 


-8 


o 'V/in'wi ron'7ono 

8.7497715807892 


10 


10.7498458 


in '7/ino/imnnn 

lU. 7498450900 


6.7 


X 


10" 


-9 


in '7/ino/irnoni/''n 

10.749845089160 


1/ 


10 '7/inooor 

12.7498885 


1 o '7 /I no o o nn o/^ 

12.7498880086 


2.0 


X 


10" 


-9 


10 '7/inooonno/ii/^ 

12.749888008416 


14 


14.7499156 


1 A '7/innicro/i/icro 

14.74991524453 


6.9 


X 


10" 


-10 


1/1 '~7 A C\C\-\ V C\ A A A A n 

14.749915244446 


lb 


1/^ '7/innooo 

16.7499oo8 


A p '7/innoo/^iinn 

16.74993361109 


2.7 


X 


1 n- 

10 


-10 


1/^ '7/innooi^iin[X'7 

16.749933611057 


18 


1 O ''7 A (\(\ A n 70 

18.74994673 


10 '7/inn/i/^ro/inr 

18.74994658405 


1.2 


X 


10" 


-10 


10 '7/inn/i/^ro/ino/i 

18.749946584034 


9n 
zu 


zu. ( ^tyyooiy 


zu. ( ^tyyoouooioi 


5.4 


X 


10" 


-11 


zu. ( ^tyyoouooi ( o 


22 


22.74996334 


22.749963259765 


2.7 


X 


10" 


-11 


22.749963259761 


24 


24.74996886 


24.749968804885 


1.4 


X 


10" 


-11 


24.749968804883 


26 


26.74997323 


26.749973181147 


7.5 


X 


10" 


-12 


26.749973181145 


28 


28.74997673 


28.749976695732 


4.2 


X 


10- 


-12 


28.749976695731 



Table 1: Numerical values of the even eigenvalues 2A„ from the large- A expansion. The first 
column gives simply the numerical values of fll.6l) . with all higher corrections ignored. 2An^'' 
are obtained from fl4.27p with the sum truncated beyond the term oc A"''. The differences 
25An'^'' = 2Ai''^ — 2An^'' are given in the third column, they show the effect of the term oc A~^. 
In the last column we present the eigenvalues 2An™™'' computed by direct numerical solution 
of (HI]). 



n 


2A„ from Eq. 01.61) 


2Ai 


2(5Ai'^ 


2_)^Cnum) 


1 


1.75381 


1.75396 


-9 


3 


X 


10" 


-4 


1.7537313369175 


3 


3.751045 


3.7510570 


-8 


6 


X 


10" 


-6 


3.7510575817054 


5 


5.750487 


5.75049257 


-5 





X 


10" 


-7 


5.7504926236487 


7 


7.7502819 


7.750284389 


-6 


4 


X 


10" 


-8 


7.7502843971925 


9 


9.7501838 


9.750185133 


-1 


3 


X 


10" 


-8 


9.7501851352539 


11 


11.7501294 


11.7501301421 


-3 


4 


X 


10" 


-9 


11.750130142515 


13 


13.7500960 


13.7500965038 


-1 


1 


X 


10" 


-9 


13.750096503972 


15 


15.7500741 


15.75007442838 


-4 





X 


10" 


-10 


15.750074428438 


17 


17.7500589 


17.75005915901 


-1 


6 


X 


10" 


-10 


17.750059159035 


19 


19.75004800 


19.750048157159 


-7 


2 


X 


10" 


-11 


19.750048157169 


21 


21.75003985 


21.750039967124 


-3 


4 


X 


10" 


-11 


21.750039967130 


23 


23.75003362 


23.750033705315 


-1 


7 


X 


10" 


-11 


23.750033705317 


25 


25.75002874 


25.750028810058 


-9 





X 


10" 


-12 


25.750028810060 


27 


27.75002486 


27.750024910393 


-4 


9 


X 


10" 


-12 


27.750024910394 


29 


29.75002171 


29.750021753287 


-2 


7 


X 


10" 


-12 


29.750021753287 



Table 2: The same as in Table 1, but for the odd eigenvalues 2A, 
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s 






(num) 






(num) 


1 

1 


2.22417142752923 


0.2241714275292 


2 


i 


^.4143983221171 


2.0000000000000 


6 


20.498 


1207536828 


1.7198 


241782619 


A 


54.538 


349992708 


1.7952480377615 


5 


147.32680373214 


1.9789 


889429098 


6 


399.32397653715 


2.2250507748184 


7 


1083.2464075913 


2.521777906136 


8 


2939.1433918727 


2.8678 


85373267 



Table 3: Numerical values of of the spectral sums (11.71) . fl4.18p . 



Again, it is easy to check that these numbers comply with (14.201) to twelve digits. 

The above numerics concern the eigenvalues {A„}. But it is also interesting to see how the 
large-A expansions in Section 4 approximate the associated eigenfunctions \Ef„(z/). It turns 
out that (14. 9 p provides a rather good approximation even if one retains only the leading 
term 1 in the expansion (14.100 of the coefficients R±{z, A). In this approximation the sum in 
(14. 7p is understood as the hypergeometric function (— i^;)^"*"^ [/(i -|- 2, — i^). This results 
in the following approximate expression for the normalized eigenfunctions (which we write 
for z/ > 0; the u < part is restored by symmetry), 



v^7rA„ sinh2(Hi) 



cosh(^)VrTe^ 
Here the phase $(z^) has expression 



•n f/(l + f , 2, -2i7rA„tanh (f ) ) 



• (5.2) 



sgn(z/) 



TT 



1 



8 An 



-2n\u\ 



2,1) 



(5.3) 



in terms of the Lerch transcendent $(z, s,a) = Yll 



The approximation is not 
u log u 



^k=0 (fc+a)= 

expected to be very accurate at small z/, because (15. 2p has a term with singularity 
at = 0, while true eigenfunctions \E'n('^) are analytic at all real u (recall that the higher order 
terms in (I4.10p were designed precisely to fix this analytic deficiency). However, numerically 
the deviations of (15. 2p from true eigenfunctions are rather small even at small u. Figs.[3lll] 
show plots of (I5.2p for low n against the corresponding eigenfunctions obtained by numerical 
solution of (12. 2p . Deviation at small u is barely visible only for n = 0. 



6 Remarks 

As was mentioned in the introduction, the techniques developed here extend to the case (II. 9p . 
While we plan to treat this case in a separate paper, let us announce here some preliminary 
results. The large-A expansion generalizes in an almost straightforward way, yielding the 
asymptotic large-n expansion of A„(a). The large-n behavior follows from the "quantization 
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Figure 4: Comparison of the approximation (15.21) for for odd eigenfuctions \l/„(z/) with 
n = 1, 3, 5 (sohd hnes) with the results of direct numerical solution of fl2.2p (bullets). 



condition", generalizing (14.271) . (I4.28P in Section 4, 



2ry ry^ 1 

2A log(2A) - Coia) + — + 



+ (I + a) 



+ 



127r8A3 
+ 0{X-Hog\\)) 



TT^A 27r6 A2 

+ vr' (1 + af - (-1)" Un^ (1 + a) (log (27re^^A) - Ci(«)) 
n, 



where 



Ci(«) 



^ + ^ log (47re^-) - 



4 7r2 

1 3a a 

2 ^ T ^ 8 



dt 



/"^ sinh(t) ( sinh(2t) - 2t) 

27r2 t cosh2(t) ( a sinh(t) + t cosh(t) ) 

sinh(2t) - 2t 
t sinh(t) ( a sinh(t) + t cosh(t) ) 

19 



(6.1) 



(6.2) 



The first two terms in fl6.ip are known since [1]. Explicit expression for tlie constant term 
Co(a), Eg. (16.21) . was previously obtained in [11] (see also [12]). We believe the higher order 
terms in the expansion in (16. ip are new. Further terms can be derived in a systematic way. 
Another result compact enough to be presented here is the analytic expression of the spectral 
sums ( 14:T# I. 



^m, X , . a f . sinh(t) (sinh(2t) ± 2t) , , 

Cl Va =log Stt -2±1-- / dt „ ;/ ^ ^^TT • 6.3 

^ ^ 4 t cosh=^(t)(asinh(t)+t cosh(t)) ^ ^ 

These and other results indicate the rich analytic structure of Xn{a) as the functions 
of complex a. First, as expected, A„(a;) have square-root branching point a = —1, which 
corresponds to the limit mi = 1712 = 0, where the chiral symmetry becomes exact. In 
particular, the lowest eigenvalue Ao(a) turns to zero as y/a + 1. But in addition, there are 
infinitely many similar square-root branching points located on the second sheet of the a- 
plane (i.e. in the left half plane of the variable y/a + 1), accumulating towards a = 00. 
At each of these points one of the even eigenvalues A2m(«) turns to zero. It is difficult to 
imagine that if one takes QCD2 with large but finite Nc these singularities just disappear. It 
is more likely that they become nontrivial critical points of some sort. What are the physics 
of these critical points? Can one identify associated (nonunitary) CFT? These are some of 
intriguing questions which we plan to study in the future. 
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7 Appendix A 

Al. Analytic expressions for the spectral sums (I1.7P with s = 2, 3, 4 are given in (II. 8p . 

(s) 

Here we present few more expressions for G'^\ with s up to 8: 

= [ ^^""^ " ^^""'^^ + ^^^^^^^ + + ^^^CiS) + 31C(5)) ] 

= [7r^(4090 + 1449C(3)) - 207r^ (l40 + 9912C(3) + 2646C(3)2 + 837C(5)) 
-M5 (138768C^(3) + 24696(^(3) + 16120C(5) + 28C(3)(140 + 837C(5)) + 3429C(7)) ] 

= — [ - 122887r6 + 497r^(21002 + 45969C(3)) - 294071^(36 + 15400C(3) 
567000 

+22050C^(3) + 4495C(5)) + 315(1509200(^(3) + 1440600(^(3) + 60760C(5) 
+196C(3)(36 + 4495C(5)) + 51689C(7)) ] (7.1) 

^ This expression follows in rather straightforward way from analysis in Appendix B. 
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(8) 



[ - 47r^(97286 + 12375C(3)) + vr^ (4210624 + 38551128C(3) 



G 

1944000 

+5622750(^(3) + 767250C(5)) - 3071^(2464 + 32104800(^(3) + 3704400(^(3) 
+4032480C(5) + 8400C(3)(712 + 279C(5)) + 142875C(7)) + 315(21403200(^(3) 
+1852200(^(3) + 76880C(5) + 432450(^(5) + 8400(^(3) (712 + 279C(5)) 
+167132C(7) + C(3)(4928 + 8064960C(5) + 285750C(7)) + 27375C(9)) ] 



G 



(5) 



G 



225 
1 



[ - 3277^ + 5071^(6 + 7C(3)) - 15(76 + 155C(5)) ] 



= — ^ [681120 - 2148007r^ + 144267r^ + 2016007r\(3) - 2 1 73 57r^C(3) 
97200 L sv ; sv y 

-1339200C(5) + 2511007r=^C(5) - 771525C(7) ] 



(7) 



G 



(8) 



G"-'' = ^ \ - 11541600 + 42453607r2 - 5193027r^ + 1843271*^ 

1190700 L 

+216097r^C(3) + 136710071^^(5) - 10921365C(7) ] 

[ 1021799520 - 42952224071^ + 603934807r^ - 2819800 



(7.2) 



76204800 

+11030880071^^(3) - 342231687r^C(3) + 14553007r*'C(3) + 259308007r^C^(3) 
-732765600C(5) + 34633200071^^(5) - 225571507r^C(5) - 3445092007r2C(3)C(5) 
+ 1144262700(^(5) - 860267520C(7) + 1260157507r2C(7) - 253519875C(9) ] 

Expressions for G^^^ with even higher s (we have them all the way up to s = 13) have similar 
structure, but appear too cumbersome to fit in reasonable page space. 

A2. Coefficients $i^^(/) in Eqs. lKTT} . Km for k<7: 



(2), 



(3), 



(4), 



(5), 



(6), 



27r4 
12/ -7 



127r6 
1 



167rS 
1 

127rio 
1 

144^ 



167r2 - 5 + 44 / - 24 /2 

3 + 767r2 + 12C(3) + (48 - 60%^) / - 84 + 24 
111 + 296571^ - 8287r^ + 948C(3) + (396 - 62287r2 - 720C(3) ) / 



+ (21607r^ - 2448) + 1968/^ - 360/ 
1 



(7.3) 



9607r 



8 70 + 40 8 6577^ - 792647r^ + 21560C(3) - 168 007r2C(3) + 4320C(5) 
- (1800 + 18168071^ - 470407r^ + 42720C(3)) / + (1603207?^ - 24240 
+ 14400C(3) ) /2 + (45760 - 336007r2) /^ - 22080 /^ + 2880 /^ 
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(2) 



(3) 



(4) 



1 

5-12/ 



127r6 
1 



1 

127rio L 
1 

144^ 



3- IGtt^ - 36/ + 24/2 
- 3 - TOtt^ - 12C(3) + (eOvr^ - 36) / + 72 - 24 
8287r^ - 87 - 26337r2 - 900C(3) + ( 58207r2 - 252 + 720C(3) ) / 



+ (1944 - 21607r^)r- 1728^ + 360/ 
1 



(7.4) 



9607r 



14 



768647r^ - 630 - 353557r^ - 19800C(3) + 168 007r2C(3) - 4320C(5) 
(1800 + 16392071^ - 470407r^ + 40800C(3) ) / + ( 18000 - 15120071^ 
14400C(3) ) + (336007r2 - 37440) + 19680 - 2880 



8 Appendix B 



Here we describe some technical details of our analysis of the integral equation (12. 8p with 
the r.h.s. (12. lip . To make the equations shorter, throughout this appendix we trade the 
variable u for 



_ nu 



but, with some abuse of notations, retain the same symbols for basic functions. Thus ^'±(t|A) 
will stand for solutions of the integral equations 



fit) 



t 



sinh(t) 

TT 

2 sinh(t) 



/oo 
dt' S{t-t') ^+{t'\X) 
-oo 

dt' S{t-t') *-(t'|A) , 



(8.2) 



oo 
oo 



with the kernel 



Sit) 



sinh(t) 



(8.3) 



The analysis below does not depend on a specific form of the function f{t). With f{t) = 
t coth(t), Eqs. (l8.2p are equivalent to (12.80 . (12. lip , but almost all statements below remain 
valid if one takes the more general form 



f{t) = a + t coth(t) 



U) 



which appears in analysis of (ll.ip with nonzero but equal ai = a2 = a. 
Eq. (l8.2p defines the spectral problem 

K(f){t) = X-^ (f){t) 



.5) 
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for the Fredholm operator 

/oo 
dt' K{t,t') (l){t') (8.6) 
-oo 

with the kernel 

where = ^/f Let R{t, t'\X) be the corresponding resolvent, i.e. the kernel of the operator 



By definition, it satisfies the equation 

/oo 
dr K{t, t) R{t, t'\X) = K{t, t') . (8.8) 
-oo 

The spectral sums (II. 7p and fl4.18p are related to the resolvent by the trace identities 

+ >^''' = C+ / dt [R{t,t\\) - R^''\t)] (8.9) 

s=l "'-^ 

°° POO 

~G^'^] X'-^ = / dtRit,-t\X). (8.10) 

-1 J —oo 



s=l 



The constant C in (18. 9p depends on the choice of the subtraction term R'^'^^t) needed to 
make the integral convergent. We take 

^(o)(,) ^ . (8.11) 

With this choice the constant can be shown to be exactly 

C = 2 log(87r) - 4 . (8.12) 

It is the remarkable property of the kernel (18. 3p in (18. 2p that the resolvent can be expressed 
in a simple way through the functions \E'+(t|A) and \E'_(t|A), namelj{§ 

i?(t,tiA) = ^/mm \^4t'm-m-^-{t'm^m] . (s.is) 

TT smh(r — t) I J 



To prove this identity, consider the Liouville - Neumann series for ^-|-(t|A), 

k 

■ atj 



°° poo -I ^ Jj. 



^.14) 



^In other words, the kernel (|8.7p belongs to the class of "integrable" kernels, see Ref. [15] for other kernels 
with similar property. 
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where to = t. Then we have 



sinh(t'^ — tr^ 



k,m 



sinh(t'^,) smh(t^ 



(8.15) 



where again to = t and t'o = t' . Let us introduce uniform notations for the integration 
variables 



[tli ■ ■ ■ tm'i tf^i ■ ■ - t^) — (Ti, . . . Tm, Tm+1, • • • Tfc+m) • 

The elementary identity 



(8.16) 



E' _sinh(r„+i - r, 



m=0 



_ sinh(r;+i - tq) 
sinh(rm+i) sinh(rm) sin(r/+i) sinh(ro) 



(8.17) 



allows one to put fl8.15p in compact form 
2 sinh(t) sinh(t' 



TT sinh(t' — t) 



fit) fit') [^4t'\x)^4t\x) -^4t'\x)^4t\x)] 



1=1 



I , i+i 



where now tq = t, r;+i = t'. It is easy to see that the right-hand side here divided by 
\/ fit) fit') is exactly the Liouville - Neumann series for the solution of the integral equation 



Now, since 



o \ p oo 

Z^±(A)=(^) exp -J2s-'G^^^y 



(8.19) 



combining Eqs.f l8.9p . (18.101) and f l8.13p leads to the following expressions for the logarithmic 
derivatives of the spectral determinants, 



dxlogiD+D^) = 2 
f ^ 



dt 



I 2/(t) 



(9a log 



Q.it\X)dtQ+it\X) - Q+it\X)dtQ4t\X) 

dt 7rQ4t\X)Q4t\X) 
^ fit) sinh(2t) 



tanh(t) 



.20) 



where 



Q±(t|A) = f sinh(t) fit) vl/±(t|A) . 



(8.211 
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The above analysis, in particular Eqs. fl8.20p . applies to (18.21) with generic f{t). If one 
takes f{t) of the special form t coth(t), it is very likely that (18.201) further reduce to the simple 
form (12.171) . Note that (I2.17P corresponds to replacing the integrals in (I8.20p . (18.200 by one 
half of the residues of the integrands at the pole at t = Unfortunately, so far we could not 
find a way to reduce the integrals to the residues, and thus (I2.17P lacks rigorous proof. But 

(s) 

it passes a number of nontrivial tests, both analytic and numerical. Thus, all listed in 
(II. 8p come out identical by direct evaluation of the integrals from (I8.20p . For higher s, using 
(12.171) instead of (I8.20p dramatically simplifies calculations, and all analytic expressions for 
G^f' listed in Appendix A and beyond in fact depend on the validity of (12.170 . We take 
agreement with the numerical data in Table 3 as further support of (12.171) . On the other 
hand, although in deriving the large-A expansion of the spectrum in Section 4 we have used 
(12.171) . it is possible to show that the results for the coefficients ^±\l) in (14.271) . (14.281) are 
independent of the validity of this relation. In particular, all expressions for these coefficients 
in Appendix A can be re-derived by a different (somewhat more complicated) method which 
does not rely on (12.171) . Let us stress also that the simplification (12.171) depends on the 
special choice f{t) = t coth(t) in (18.21) . It is unlikely that any simple modification of (12.171) 
exists for more general f{t), say of the form (18.41) . Therefore, in the analysis of the problem 
(II. ip in the more interesting case of a generic a (which we plan to present in a separate 
paper), we have to make do with the integral representation (I8.20p . 
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